Spinmotive Force due to Intrinsic Energy of Ferromagnetic Nanowires 
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We study, both analytically and numerically, a spinmotive force arising from inherent magnetic 
energy of a domain wall in a wedged ferromagnetic nanowire. In a spatially-nonuniform nanowire, 
domain walls are subjected to an effective magnetic field, resulting in spontaneous motion of the 
walls. The spinmotive force mechanism converts the ferromagnetic exchange and demagnetizing 
energy of the nanowire into the electrical energy of the conduction electrons through the domain 
wall motion. The calculations show that this spinmotive force can be several microvolts, which is 
easily detectable by experiments. 



Recently, an electromotive force of spin origin has been 
theoretically predicted [1[ and experimentally observed [2|, 
y] . This "spinmotive force" reflects the energy conversion 
of magnetic energy of a ferromagnet into electric energy 
of conduction electrons via the spin exchange interaction. 
This offers a new functionality for future spintronic de- 
vices. Theoretically, it has been pointed out that the 
spinmotive force is produced by the spin electric field 
acting on conduction electrons [j4S| • 
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where m is the unit vector paralell to the direction of the 
local magnetization in a ferromagnet, P the spin polar- 
ization of the conduction electrons, and e the elementaly 
charge. It is requied that the magnetization depends 
both on time and space. To demonstrate the spinmo- 
tive force, Yang et al. studied a field-induced domain 
wall (DW) motion in a uniform nanowire [2], in which the 
moving DW releases its Zeeman energy into conduction 
electrons through the spinmotive force mechanism. On 
the other hand, it has been suggested that a DW motion 
can be induced solely by "shape effect" [9, 10]. In addi- 
tion to the Zeeman energy, a DW has a surface-tension- 
like intrinsic energy due to the ferromagnetic exchange 
and magnetic anisotropic energy of the magnetization. 
Thus in a spatially-nonuniform nanowire a DW sponta- 
neously tends to move to narrower regions in order to 
reduce its intrinsic energy, which is proportinal to the 
cross-sectional area of the wirejol, |lO|. The shape-effect- 
induced DW motion thereby gives rise to a spinmotive 
force without any external field, using the inherent energy 
of the magnetic texture, i.e., the ferromagnetic exchange 
and anisotropic energy [9]. In reality, the DW motion 
in a spatially-nonuniform nanowire shows highly nonlin- 
ear behaviors |lCHl4j . for a quantitative understanding of 
which numerical methods are of particular importance. 

In this paper, we study, both analytically and nu- 
merically, the spinmotive force in a wedged magnetic 
nanowire. The results of the numerical simulations 



strongly support this concept, and show that this spin- 
motive force can be large and lasts long enough to be 
detected in experiments. 

We consider a Permalloy thin wire in which one of 
its ends is wedged as illustrated in Fig. [H The length, 
width, and thickness of the wire are 5 /im, w = 100 nm, 
and d = 4 nm, respectively. The cross-sectional area of 
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FIG. 1: (a) The schematic illustration of a NisiFeig wedged 
nanowire containng a single DW. The arrows indicate the 
magnet izaiton direction m. The DW is characterized by two 
collective coordinates, i.e., q, the center position, and ip, the 
tilt angle from the easy plane, (b) The spinmotive force and 
the shape-effect magnetic field, calculated by Eqs. (|2]), Q, 
(O, and Q, with if = and the tail-to-tail configuration, 
as a function of the DW position q. The inset describes the 
dimensions of the sample. 



the wire A tapers toward the end as 



A{z) = ( --z^w ] d. 



(2) 



First, we analyticaly investigate influence of the wedge 
on the spinmotive force based on a cohective-coordinate 
model. Here we employ the one-dimensional magnetiza- 
tion dynamics along the z-axis, assuming uniformity in 
the lateral direction (the x-y plane), in which the dynam- 
ics of a DW is characterized by two collective coordinates, 
i.e., the center position q and the tilt angle t/j of the DW 
plane relative to the x-z plane as indicated in Fig. [T](a). 
The time evolution of {q^ tjj) is described by the reduced 
Landau-Lifshiz- Gilbert equations [101] 
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where 7 is the gyromagnetic ratio, a is the Gilbert 
damping constant, A is the wall width parameter, H 
is an external magnetic field, and Hk is a perpendicu- 
lar anisotropy field due to the demagnetizing interaction. 
Here and hereafter, upper (lower) sign corresponds to the 
tail-to-tail (head-to-head) configuration of the DW. In a 
spatially- nonuniform nanowire, in addition to H^ DWs 
are subjected to the "shape-effect" field [ol, llO|: 



H. 



shape 



= ± 



d 



2Msdq 



InAiq), 



(5) 



where a is surface tension energy of a DW and Mg is 
the saturation magnetization. As in the text book[15j, 

energy minimization gives a = ctq (l + Q~^ sin^ 7/^) 
Here ctq = \^^A^K^ with A^ the exchange stiffness con- 
stant, Ku the uniaxial anisotropy constant, and Q is 
the ratio of K^ to the perpendicular anisotropy con- 
stant. This illustrates that the gradient of the intrinsic 
DW energy (jA{q) exerts a force acting on the DW. In 
the present system, the magnetic anisotropy comes from 
the demagnetizing field while the crystalline magnetic 
anisotropy is negligible. It is the exchange and the de- 
magnetizing interaction between the magnetization that 
is responsible for i^shape- 

From Eqs. ([3j) and (j4|), in a spatially- nonuniform mag- 
netic nanowire, i^shape can drive a DW and make precess 
the DW plane even without the external field H. In the 
light of this simplification, with some sign convention the 
spinmotive force across the DW is obtained as 
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V= I dzE, = ± — V^. 
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Equation (J6]) indicates that the sign of the spinmotive 
force reflects that of the precession rate of the DW plane 
-0, and is independent of the translational motion of 



the DW[ll, [3. When | i:^ + i:^shape | is smaller than the 
Walker breakdown field i^w, the Larmor torque due to 
H + i^shape cannot overcome the countertorque due to 
the perpendicular anisotropy; H + i^shape in Eq. (|4]) 
are canceled out by the third term with a certain i/j, re- 
sulting in ip = 0. In contrast, if \H -\- H chapel is larger 
than i^w = (^Hk/2^ the DW plane precesses with fi- 
nite ip. Thus, the spinmotive force appears only when 
\H + i^shapel is abovc the Walker breakdown field. Once 
the DW configuration is fixed, the sign of the spinmotive 
force is determined by the direction of H -\- i^shape- 

The shape-effect field i^shape and the spinmotive force 
V, calculated from Eqs. (|2|), ©, dSJ and (O, with H = 
and the tail-to-tail configuration, as a function of the 
DW position q are shown in Fig. [1] (b). Typical values 
for a Permalloy nanowire are employed; 7 = 1.76 x 10^^ 
Hz/T, a = 0.01, Ms = 1 T, As = 1.3 x 10"^^ J/m, 
K^ = 10^ J/m^, and P = 0.6. We neglect the third 
term in Eq. (JD), in which aHx is tens of Oersted for 
the Permalloy nanowire, and a c^ ctq is assumed. In this 
scheme V is proportional to i^shape, and in the wedged 
part {z > 0) the finite spinmotive force appears. Figure [1] 
(b) indicates that the spinmotive force of several micro- 
volts can be obtained in the present system without any 
applied field but using the intrinsic magnetic energy of 
the nanowire. 

The collective-coordinate model imposes a number of 
assumptions, i.e., the uniformity in the x-y plane , a rigid 
wall structure, and constant magnetic shape anisotropics. 
In fact, these assumptions are no longer guaranteed with- 
out check for nanowires with the shape modulation and 
in particular for relatively higher driving field region 
{\H -\- Hshapel :^ Hw)- To confirm the above predic- 
tion and to obtain the detailed real-time spectrum of the 
induced spinmotive force with high accuracy, here we per- 
form numerical calculations free from the constraint. We 
solve the original Landau-Lifshiz-Gilbert equation 

rh (r , t) = —7m (r , t) x ffeff + <^^^ (^, t) xrh (r , t) , (7) 

using the OOMMF code [16] with dividing the sample into 
4x4x4 nm^ cells. Here ffeff is the effective magnetic 
field including the external, exchange and demagnetizing 
fields. The parameters are the same as used in the one- 
dimensional calculation. Here the demagnetizing field is 
taken into account by the dipole interaction, instead of 
K^. To track the DW position we calculate the time evo- 
lution of spatially- averaged value of the 2:-component of 
the magnetization, (rriz) = {1/V) J rrizdV^ as shown in 
Fig. [2] (a). Initially (t < 0), a transverse tail-to-tail DW 
is trapped at z = — 1 /im, where the 4x4x4 nm^ defect 
is introduced. From t = to 1.9 ns a constant magnetic 
field pulse H = (0, 0, H), where H = —20 Oe, is applied, 
see Fig. [2](b); during this period, the DW has been dis- 
lodged from the defect, traveled about 1 /im and entered 
the wedged part. The DW maintains its transverse struc- 
ture through the field-induced motion. Even after the 
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FIG. 2: The time evolution of (a) the spatially averaged z- 
component of the magnetization {rriz) (the inset showing a 
magnified image) with the five snapshots of the distribution 
of ruz corresponding to time indicated by the arrows, (b) the 
applied magnetic field, and (c) the spinmotive force between 
the ends of the nanowire. 



applied field is turned off at t = 1.9 ns, the DW keeps 
moving forward and sHdes down the slope with the help 
of the shape-effect. In this stage there arises an oscilla- 
tion of the DW position, i.e. the back and forth motion of 
the DW, with the periodic nucleation and annihilation of 
ant ivort ices. This is because of the Walker breakdown in- 
duced by the shape-effect. In the one-dimensional model, 
as the tilt angle '^ develops through the Walker break- 
down motion, (A7J^x/2) sin2'0 in Eq. ([3j) oscillates and 
thus q can be negative. This oscillatory motion contin- 
ues until the DW reaches the end of the nanowire and 
disappears at t ::^ 20 ns. We remark here that after the 
whole event, the exchange and the demagnetizing energy 
of the nanowire has obviously decreased. The spinmotive 
force mechanism, which guarantees the total energy con- 
servation, dictates that the decreased magnetic energy 
is received by the conduction electrons, resulting in an 
electromotive force. 

A measurable voltage between given two points can 



be calculated as the spatial difference of the spin- 
independent electric potential ^(r,t) between the points, 
which is obtained by sol ving the Poisson equation 
V^(/)(r,t) = -V • E^{r,t)\^,^. Figure [2] (c) shows the 
spinmotive force between the ends of the nanowire [see 
Fig. [U (a)] versus time. In spite of the field-induced DW 
motion, from t = to 1.9 ns, a very small spinmotive 
force is observed since |JY| = 20 Oe is below the Walker 
breakdown field for the present nanowire. After the sys- 
tem goes into the Walker breakdown, with no applied 
magnetic field, the finite spinmotive force is induced un- 
til the DW vanishes at t ::^ 20 ns. The observed spin- 
motive force is of the order of microvolts, oscillating but 
unipolar. These facts agree with the analytical predic- 
tion from Eq. (J6]) . The oscillation of the spinmotive force 
reflects the DW structure transformation. Apparently il) 
becomes larger in an anti- vortex motion than in that of 
a transe verse wall and so does V . The peaks of the spin- 
motive force signal in Fig. [2] (c) are associated with the 
traverse motion of an antivortex, i.e., the plateau in the 
DW position curve in Fig. [2] (a). The DW moves back 
and forth with changing its structure, whereas, the sign of 
the spinmotive force is unchanged and determined by the 
direction of the shape-effect magnetic field as estimated 
in Fig. [T](b) [see also the discussion below Eq. (|6])]. 

In conclusion, we have investigated the spinmotive 
force in a wedged magnetic nanowire, and numerically 
confirmed the spinmotive force of several microvolts for 
tens of nano seconds, using only the inherent magnetic 
energy of the nanowire. The behavior of the spinmotive 
force depends on the shape and the material nature of 
the sample, which implies the posibility of various appli- 
cations. 
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